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Abstract 



We consider two parameter families of differential equations on a Banach space X, 
where the parameters c and e are such that: 

• when e = 0, the differential equations are symmetric under the action of the 
group of one-dimensional translations SE(1) acting on X, whereas when e 7^ 0, 
this translation symmetry is broken, 

• when e = 0, the symmetric differential equations admit a smooth family of rel- 
ative equilibria (travelling waves) parametrized by the drift speed c, with c = 
corresponding to steady-states. 

Under certain hypotheses on the differential equations and on the Banach space X, we 
use the center manifold theorem of Sandstede, Scheel and Wulff [21] to study the effects 
of the symmetry-breaking perturbation on the above family of relative equilibria. In 
particular, we show that the phenomenon commonly referred to as propagation failure, 
or wave blocking occurs in a cone in the (c, e) parameter space which emanates from 
the point (c, e) = (0,0). 

We also discuss how our methods can be adapted to perturbations of parameter- 
independent differential equations (such as the Fisher-KPP) which admit families of 
relative equilibria parametrized by drift speed. 
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1 Introduction 



Travelling waves are an important class of solutions for many types of systems which are 
modelled using reaction-diffusion partial differential equations [31 EE [91 [121 E31 EH EH 123] . 
integro-differential equations [21 [151 E2] or delayed partial differential equations [31 El [3 EHj . 
In biological systems, travelling waves frequently represent electrical pulses, fronts, backs or 
wave trains which propagate along one-dimensional networks of excitable cells [H1[T8]. They 
have also been observed in chemical reactions, combustion theory [23] and nonlinear optics 
PQ. In fact, the literature on travelling waves, both theoretical and applications, is quite 
extensive. We will however single out the excellent review articles by Sandstede [2Q] and by 
Xin [25]. 

For nonlinear models of propagation in spatially extended systems, a common modelling 
hypothesis or simplifying assumption is to assume that the underlying medium of propaga- 
tion of the wave is homogeneous. However, this hypothesis is seldom reasonable in many 
practical applications. Indeed, biological media typically contain many sources of inhomo- 
geneity (e.g. eschemic tissue, gap junctions), which may have effects on the propagation of 
the wave through the medium (i.e. acceleration, deceleration, reflection, or blocking propa- 
gation) [31 [T21 EEZ1 1231 [2S]. It is this interplay between inhomogeneity of the medium and the 
propagation properties of travelling waves which we seek to elucidate in this paper, from a 
model-independent point of view. 

Mathematically speaking, travelling waves are typically a consequence of translation- 
invariance (homogeneity) of the underlying mathematical model. To illustrate this point, 
consider the following well-known bistable reaction-diffusion equation (see e.g. [H]) 

9u d 2 u ^ . . . 

m = S? + -«)(«-«) d-D 

where u(x, t) is a function of x G R and t > 0, and < a < 1 is a real parameter. Equation 
(II. ip is invariant under the change of variables 

v{x,t) =u(x + 7,t) (1.2) 

for any 7 G R. Using the language of group-equivariant dynamical systems [8], we say that 
( 11. 21) defines an action of the additive group of real numbers SE(1) on a suitable space of 
functions, and that the right-hand side of equation ( II. ip is equivariant under this action. One 
consequence of this symmetry property (see e.g. [8]) is that solutions of ( II .ip are mapped 
to other solutions by elements of the group SE(1). A travelling wave solution to (II. ip is a 
solution such that its time orbit is contained in its SE(l)-orbit. Such a solution is called a 
relative equilibrium. 
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To find relative equilibria in (II. ip . we assume the ansatz 

u{x,t) = u*(c){x + ct), (1.3) 
which reduces (II. ip into an ordinary differential equation for the function u*(c)(r]): 

du*(c) d 2 u*(c) 



drj drf 



+ u*(c)(l-u*{c)){u*{c)-a). (1.4) 



One then uses techniques from phase-plane analysis to seek values of c for which ( 11. 4p 

/ du*(c)\ 

admits a heteroclinic orbit connecting the equilibrium points I u*(c), — - — = (0,0) and 



du*(c) 

u*(c), — ) = (1,0). In particular, one can show [TJ] that 



dr) 



drj 



u*(c)(rj) — - ( 1 + tanh 



2 V \2V2 
is such a heteroclinic solution to (II. 4p with 

c=-L(l-2a). (1.5) 

Using the SE(l)-equivariance, it follows that for any 7 G R, 

u(x,t) = M*(c)(x + 7 + ct) = ~ ^1 + tanh { ^ + ^ I 1 - 6 ) 

is a travelling wave solution to ( II. ip for c given by (II. 5p . These solutions are such that 
u — y 1 as x — y 00 and u — >0asx— y — 00, and represent a wave front which propagates from 
right to left when c > (a < 1/2) and a wave back propagating from left to right when 
c < (a > 1/2). Note however that when a = 1/2, we have c = in (II .5p . and thus (II. 6p 
represents a one-parameter family of steady-state solutions 

u(x,t) = u*(0)(a; + 7) = - ( 1 + tanh ( ^2 ) ) V7 E R. 



2 V V 2^ 

The above analysis of (11.11) is standard for many classes of equations which admit travelling 
wave solutions in homogeneous media, i.e. suppose that the travelling-wave ansatz (II. 3p 
holds, and seek homoclinic (for pulses) or heteroclinic (for fronts and backs) solutions in a 
reduced equation. 
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Suppose now that we perturb fll.ip by introducing a small non-homogeneous term 



~q[ = q^2 + u ( 1 -u){u-a) +eg(x), (1.7) 

where g(x) is some smooth (non-constant) bounded function, and < e C 1 is a small 
parameter. Equation ( 1 1.7ft no longer admits the translation symmetry of f ll.ip . so pure trav- 
elling wave solutions such as (JL3]) are unlikely to exist. However, by continuous dependence 
of the solution on parameters, we expect solutions of fll.7j) which are "close" (in some sense) 
to (jl. 6p for e small enough. There are essentially two regimes in parameter space to consider: 
Regime I. If the computed value of c from formula (II. 5p is non-zero and large relative to e, 
then heuristically speaking we expect to have a solution of ( 11.71) which is approximately a 
wave front whose "propagation speed" is spatially-dependent (c(x) ~ c), but non-vanishing 
on BL In other words, in this regime, we expect travelling wave solutions of ( 11.11) to persist 
as solutions of (II. 7p with wave speed modulated by the heterogeneity, but the perturbation 
is not large enough to stop the wave from propagating along the real line. 
Regime II. When c and e are of the same order and are both small, propagation of the wave 
through the inhomogeneity may no longer be possible because of the possible persistence 
of steady states of ( 11. ip which may survive the perturbation in ( 11. 7ft . This phenomenon 
is known as propagation failure, or wave blocking, and has been extensively studied in the 
literature for various specific models (e.g. [3j [121 HHJ EH [231 125]). 

Our point of view is that most of the above discussion which was centered around the 
specific example ( II. ip and ( II. 7p does not depend qualitatively on the particular differential 
equation, but only on the fact that translation symmetry is broken by the addition of a small 
perturbation. We will thus consider a family of abstract differential equations parametrized 
by two parameters: c and e. The parameter e controls the translational symmetry-breaking, 
while the parameter c parametrizes the drift speed of relative equilibria (when e = 0). Our 
main result will be to show that the point (c, e) = (0, 0) in parameter space acts as an 
organizing center for the phenomenon of wave-blocking. Indeed, we will show that wave- 
blocking occurs in a cone emanating from this organizing center, and that this is a universal 
feature for any abstract differential equation which undergoes forced translational symmetry- 
breaking (under certain hypotheses). 

To our knowledge, the study of the phenomenon of wave-blocking using techniques from 
dynamical systems for reaction-diffusion equations was first undertaken in [T7]. In that 
paper, the approach was to use a combination of comparison methods, bifurcation theory 
and numerical analysis to study a fixed travelling wave (c > 0) in scalar reaction-diffusion 
equations, and let the size of the inhomogeneity increase until a saddle-node bifurcation 
occurs, and the resulting equilibria block propagation. In this sense, the analysis in [17] is 
focussed away from the above-mentioned organizing center at (c, e) = (0,0). 
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Recently, in [23], an analysis for a specific inhomogeneous system of react ion- diffusion 
equations was undertaken, using a formalism which is somewhat related to the approach 
we present here, i.e. centering the analysis around the codimension 2 drift bifurcation. 
However, our approach will be general enough to include reaction-diffusion equations, as well 
as integro-differential equations, delay-reaction-diffusion equations, and potentially more. 
In [T5], an analysis of wave-blocking in a specific system of integro-differential equations 
for inhomogeneous neural networks was undertaken, and equations for wave-blocking were 
derived which are similar in spirit to the equations we will derive here for general systems. 

Our approach will be to use the center-manifold theorem of Sandstede, Scheel and Wulff 
[2T] for abstract group-equivariant evolution equations on Banach spaces as a key ingredient 
in our analysis. Specifically, we consider the semi-linear equation 



which has SE(l)-symmetry if and only if e — (i.e. Q breaks the translation symmetry). 
When e = this equation will admit a smooth two-dimensional invariant manifold of solu- 
tions foliated by travelling waves (relative equilibria) with drift speed c G (— Cq,Cq) for some 
cq > 0. Assuming normal hyperbolicity, this manifold persists when e ^ 0, but the dynamics 
on this manifold are perturbed. In particular, the line of equilibria c = for the unperturbed 
equation is deformed, and it is this perturbed family of equilibria which causes propagation 
failure when (c, e) lies in some cone in parameter space which emanates from the point (0, 0). 
Our main result to this effect is Theorem 14.31 

The paper is organized as follows. In section 2, we give the functional analytic group- 
equivariant framework for our analysis. In section 3 we consider the unperturbed (i.e. pos- 
sessing full translation symmetry) differential equations, and apply the center manifold the- 
orem of Sandstede, Scheel and Wulff [2TJ. The curve corresponding to zero drift speed on the 
center manifold will be central to the analysis of the perturbed equations, which is done in 
section 4 using the global parametrization near relative equilibria presented in [21]. In section 
5, we briefly discuss modification of our method to study parameter-independent equations 
which admit families of travelling wave solutions such as the Fisher-KPP equation. We end 
with some concluding remarks in section 6. 



2 Functional analytic setup and translation symmetry 





Let X be a Banach space with norm 
equation on X: 




and consider the following nonlinear differential 
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where c and e are real parameters, e small, A : X — y X is a sectorial densely defined closed 
linear operator, the functions J 7 : X a x R — y X and Q : X a x R — y X are C k+2 smooth 
for some k > 1 and some a £ [0, 1) (see [H]), and is uniformly bounded on X x [0, Eq\ for 
some £o > 0. Under these conditions, (12. ip generates a local semi-flow $t(u; e, c) on y := X a 

m 

Now, let T = SE(1) denote the group of rigid translations of the real line, which is 
isomorphic to the additive group of real numbers. Let T : T — y GL(Y) be a faithful 
strongly continuous and isometric representation of T in the space of bounded invertible 
operators on Y. For a £ T, we will use the notation T a to denote T(a). 

The infinitesimal generator of T is the densely defined linear operator £, defined by 

£u = lim -(T a - I)u, 

wherever this limit exists. 

Hypothesis 2.1 We make the following assumptions on the functions J 7 and Q and the 
operators £ and A: 

(a) The domains of the linear operators A and^, T>(A) = Y\ andV(C, ) = Yi, are both dense 
and V -invariant. Furthermore, there is a dense V -invariant subspace Y3 C liflY^ such 
that A(Y 3 ) C Y 2 , £(y 3 ) C Yi, and£A = A£ on Y 3 . 

(b) AT a = T a A onYt. 

(c) The function J 7 is T-equivariant, i.e. 

T(T a u, c) = T a T(u, c), Vu £ Y, c £ R (2.2) 

This implies also 

D u T(J~ a u,c)T a v = T a D u J z (u,c)v, Wu, v, £ Y, c £ R, (2.3) 

^ c) = D u T(u, c)iu V« £ F 2 , c £ E, such that J*(w, c) £ F 2 , (2.4) 

D c T(T a u, c) = T a T> c T(u, c), Vn £ Y, c £ R (2.5) 

where DuJ 7 and denote respectively the partial Frechet derivatives of J 7 with re- 
spect to u and to c. 

(e) The function Q is not T-equivariant 
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We say that (12. ip undergoes forced translational symmetry-breaking when e 7^ 0, since 
the equation admits T as a symmetry group if and only if e = 0. 

Here are a few examples of applications where Hypotheses 12.11 are satisfied. 

Example 2.2 Consider the following nonlinear reaction- diffusion equation 

du d 2 u 

—(x,t)=D—(x,t) + f(u(x,t),c) + eg(x), (2.6) 

where u(x) = (ui, (x), . . . , u n (x)) T , D is a diagonal n x n matrix with positive entries, 
f : R™ x R — y R n is smooth, and g : R — > R n is smooth, non- constant and bounded. 

We set X = £ 2 (R,R") and Y 1 = T> ( ^— ] = H 2 (R,W l ). The action ofT = SE(1) on 

\ox 2 J 

Y = X a is defined by 

(T a u)(x) = u(x + a). 
The infinitesimal generator for this action is 

(£ u)(x) = u'(x) 

with dense Y -invariant domain 

Y 2 = £>(£) = H l (R,R n ). 

Equation Ii2.6\) is invariant under this action if and only if e = 0. For the subspace I3 of 
Hypothesis \2.1\ (a), we can choose 

Y 3 = H 3 (R,R n ) 

and can easily verify that ^.(^3) C Y%, ^(^3) C Y\ and At; = C,A on Y3. 

Example 2.3 Let A be a constant nxn matrix, and consider the following nonlinear integro- 
differential equation which arises frequently in neural field models with nonlocal interactions 

m 

du f°° f°° 

— {x,t)=Au(x,t) + j W(x-y)f(u(y,t),c)dy + e J D(y)g(u(y,t)) dy (2.7) 

where W{rf) and D{jf) are smooth nxn matrix-valued functions, D is non-constant, and 
/ : M" x R — y R" and g : R™ — y R n are smooth and bounded. Furthermore, W , D, f and g 
are such that the integrals in \2. 7| ) are all well-defined. We have X = Y = C°^/R, R n ) , and 

the SEi(l)-action is given by translation of the domain, as in Example \2.2\ above. Equation 
\2. 7| ) has T = SE(1) symmetry if and only if e = 0. We set 



Y 1 = Y, F 2 = y 3 = ^ n ./R,R"). 
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Example 2.4 Let % be the Hilbert space £ 2 (R, R) and consider the Banach space X = 
C°([— r, 0], H), r > 0, with sup norm. For a function u : R — > H, define u t G X by 
u t (6)(x) = u(t + 0)(x) for i 6 i and 9 G [— r, 0]. Suppose A : X — > H is a bounded linear 

operator, A : T>(A) C H — > H is the densely defined operator Au{t) = -^-^(t,x), where 
u(t)(x) = v(t, x). Consider the retarded partial functional differential equation 

u{t) = Au(t) + A(i*t) + F(ut, c)+e G{u t ) (2.8) 

where the nonlinearities F:IxR — y H and G : X — > H are sufficiently smooth. 

The group SE(1) acts on X using the standard action on H as in Example \2.2 above: 

T a (u(6)(x))=u(6)(x + a). 

We suppose that A and F are such that 

A(T a <p)(x) = A(<p)(x + a), V^Gl.Vael 

and 

F(T a <p,c)(x) = F(tp,c)(x + a), Vy> Gl.VaG R, VcG R, 

but G does not satisfy this symmetry property. We also assume that ifipE C°([—r, 0], 
then A(<p) G if 1 (R, R). We note that these conditions on A and F are satisfied in many ap 
plications. 

As was done in J7^/ . we define 



x _ x m _ / In if 6 = 



and consider the Banach space BX of functions from [— r, 0] to X which are uniformly 
continuous on [— r, 0) and wtt a jump discontinuity atO. Elements of BX can be written as 
U = u + Xqci with u G X and a G %, so that BX is identified with X x H. The differential 
equation A2.8\) can then be written as an abstract ODE on BX: 

uv 

— = Av + X (F(v,c)+eG(v,x)) (2.9) 
at 

where A : BX — > BX is defined by 

Atp = ^ + x [A(4j) + A4j(o) - v>(o)], 

on the domain 

V(A) = Y x = {ip G X : tp G X, <p(0) G H 2 (R, R)}. 
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The action o/SE(l) on BX is: 

T a (u(8)(x) + X a(x)) = u(9)(x + a) + X a(x + a). 

Equation A2.9\) has SE(1) symmetry if and only if e = 0. The infinitesimal generator £ of 
SE(1) is 

£ (u{9){x) + X a{x)) = u{6)\x) + X a'(x) 



defined on 



£>(£) = F 2 = C°([-r, O^-H^R)) x H\R,R) cX xH = BX, 



where the ' denotes a derivative with respect to x for an element of H l (R, M.) . We may then 
choose 

Y 3 = {ip E X : ip, up E C°([-r,0},H l (R,R)), <p(0) E H 3 (R,R)}. 

It is then easy to verify that Y\, Y 2 and Y 3 are SE(1) -invariant, and that £,(Y 3 ) C Y\, 
A(Y 3 ) C Y 2 and that £A = A£onY 3 . 



3 Relative equilibria in unperturbed equations 

For our analysis, it will be useful to single out the parameter c in (12. ip and consider a 
suspended system 



^ = Au + F(u, c) +eQ(u,e) 



dt 



(3.1) 



on the Banach space Y x R, with norm || ||y X R defined by 



\u\ 



c . 



Yxl 



We extend the action T a of V to Y x E using the trivial representation on the R factor, i.e. 



u 



TaU 



' a 
C 
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If £ = in (13. ip . then this equation is r-symmetric. The V symmetry is broken when 
e 7^ 0. In this section, we will study ( 13. ip with e — 0, i.e. 

(in . . 

— = Au + s{u,c) 
at 

(3.2) 

dc 

— = 0. 

For a given fixed value of c, a relative equilibrium for ( 13. 2p is a solution of the form 

C (t) ; - ^ c 

where 76I, and u* G F3 is such that ^(u*, c) £ Y 2 . Substitution into the first equation of 
(13. ip and using equivariance yields 

Tytliu* =% t (Au*+ JF(u*, c)). (3.3) 

If we assume that ^ 0, it follows that 

= N^ + M , (3.4) 

II 

The real number 7 is called the drift speed of the relative equilibrium. We now make the 
following hypothesis for ( 13. 2ft : 

Hypothesis 3.1 We assume that there is an interval (c a ,Q>) C R and a smooth function 

u* : (c a , c 6 ) — > Y 3 

such that 

(a) For all c G (c a , c&) , we /iai>e J 7 (it*(c), c) G Y2, £ w*(c) 7^ 0, and both ^ u* (c) and Au*(c) 
are smooth functions of c, with derivatives given respectively by £w*(c) and Au*(c). 
Furthermore, we assume that T a u*(c) = u*(c) if and only if a = (i.e. u*(c) has trivial 
isotropy). 

(b) For all c G (c a ,c&), k3. 1\) has a relative equilibrium 

u(t) \ _ f %( c ) t u*(c) 
c(t) ) \ c 

where the drift speed 7(c) is the smooth function defined by 

7(c)£u*(c) =Au*{c) + JF{u*{c),c). (3.5) 
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(c) The function 7(c) defined above is monotone on (c a ,Q,) ; and there exists Cq G (c a ,Cb) 
such that 7(c ) = 0. 

Proposition 3.2 Without loss of generality (upon changes of variables and relabelling of 
functions) we may assume that 7(c) = c in 113. 5\) . i.e. 

c£u*(c)=Au*(c)+J r (u*(c),c), VcnearO. (3.6) 

Proof The function 7(c) is invertible near c = Cq, so we may define c = 7(c), u*(c) = 
m*(7~ 1 (c)) and J-{u, c) = F{u, 7~ 1 (c)). Equation (13. 5 p then becomes 

5Zu*( 1 - 1 (5))=Au*( 1 - 1 (5))+F(u*( 1 - 1 (5)), 1 - 1 (5)), 



or 



c f m* (c) = .Aw* (c) + ^(w* (c) , c) . 
We get (I3.6p by relabelling c — >■ c, u* — » w* and J 7 — )■ J 7 . I 

In the sequel, we will assume that this change of parameter and relabelling of functions 
has already been made in (I3.2p and (13.11) . 

Setting c = in (I3.6P and applying £ yields 

A£u*(0) + D u F(u*(0),0)tu*(0) = (3.7) 

Also, differentiation of (13. 6p with respect to c and evaluating at c = yields 

Au* c (0) + D u F(u*(0) } 0) u*(0) + D c T(u*(0), 0) = £ u*(0) (3.8) 

Thus, the following densely defined closed linear operator 

L:Y xR — >Y xR 



L 

is such that 



<p \ = ( Aif + DJF(u*(0),0) ^ + D c T(u*(0),0) uj 
w " 



(3.9) 



i.e. L has an eigenvalue at with eigenvector <pi = (£ u*(0), 0) T and generalized eigenvector 
<p 2 = (u*(0), 1) T . We make the following further assumptions: 
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Hypothesis 3.3 We will assume that 

(a) The eigenvalue atO of L is such that the corresponding generalized eigenspace E is two- 
dimensional, i.e. E = span{0!,02}. Furthermore, all other elements of the spectrum 
of L are bounded away from the imaginary axis (spectral gap condition). We assume 
that the Fredholm index of L is 0, and that we have the following L-invariant splitting 

YxR = E(BW = span{0i, </> 2 } © W (3.10) 

with bounded projection P, range(P) = E, ker (P) = W , and L(W) = W. We will 
also use the projection 

Q : E — > E 

such that range(Q) = M{02}- Note that since any element of W = L(W) must have 
the form 



£u*(0) 

and since <pi = \ J > then 



n 1 = 1 n 



J 10 

(b) The map a ^ T a u*(0) is C k+2 inaeT. 

(c) For any k > 0, there exists a 5 > such that ||7^u*(0) — u*(0)|| > 5 whenever \a\ > k. 

(d) The map a i— > T a v is C k+l whenever v G E. 

(e) If P denotes the spectral projection onto E in item (a) above, then the projections 
T a PT~ a are C k+1 smooth in a G T in the operator norm. 

Proposition 3.4 The mapping u*(c) in Hypothesis \3.1\ can be chosen so that 

u*(c)-u*(0) \ _ ( 



P ( v ; y M = ( Q ) , for all c near 0. (3.11} 



13 



Proof Clearly QP ( J j = 1 ' ^ S6e ^ e en< ^ °^ ^ em ( a ) * n Hypothesis 



71 a u*(c)-u*(0) 




w*(c) -u*(0) ^ _ ( 
( 

3]). Consider the smooth map 

/ : R x I 

defined by 

f(a,c) = {I-Q)P 
We have /(0, 0) = and 

/ a (O,O) = -(/-g)P0 1 = -0 1 ^O. 

By the implicit function theorem, there exists a smooth function a(c) such that a(0) = 
and such that 

(/ - Q)P ( T ^ ^ ~ M * (0) ) = ( ° ) , for all c near 0. 

We then redefine 71 a ( c ) u*(c) — ► u*(c). I 
Remark 3.5 It follows from Proposition^^ that 



Hypothesis 13.31 allows us to use the center manifold theorem of Sandstede, Scheel and 
Wulff [21], from which it follows that the set 

S — t I *^~ aU ^ ) ; oGT.c near 



is a normally hyperbolic C fc+1 smooth invariant center manifold for ( 13. 2p . The dynamics 
of (13. 2p restricted to S is rather trivial, described by the two-dimensional center manifold 
ordinary differential equations 

t : o. < 3 - i2 > 

In particular, the line of equilibria c = for ( 13 . 1 2[) corresponds to the smooth one- dimensional 
group orbit 



e-{(^);. 6 r} 
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of the steady state I U ^ ) of (lOl . 





As was shown in [2T], we obtain a global parametrization of a neighborhood of as 
follows: any point u close enough to G can be written as 

u = U\(u) + U 2 {u) + w(u), 

where 

UM = T a{u) ( M ^ 0) ) , U 2 (u) G r a{u) (range(QP)), w{u) G T a(u) (W), 

where all functions of u above are C k+l . In particular, for any c close enough to 0, we may 
write the following representation of the center manifold S: 



r a u*(c) 

c 



(3.13) 



where by Proposition 13.41 and Remark 13.51 we have 

v{c) \ _ f u*(c) - u*(0) - cu*(0) 

x J' V 

i.e. 

v(c) \ _ ( 



G W, 



/ V 

where P is the projection operator onto E in Hypothesis 13.3( a). See figures [T] and |2] for a 
schematic representation of S and of the flow of f)3.12p on S. 

Taking into account the equations (I3.12p on the center manifold, the decomposition (13. 13[) 
when substituted into (13. ip leads to the compatibility condition 

^p(^f ') + cp(^' W)=p( 7/ ( C ) ) (3.14) 

where 

H(c) = T(u*(0) + cu* c (0) + v(c),c) - T(u*(0),0) - cD u T(u*(0),0)u* c (0) 

-cD c JR(u*(ti), 0) - D u F(u*(0), 0)v(c). 

Using normal hyperbolicity arguments similar to those of LeBlanc and Wulff [16], the 
set S will persist as a normally hyperbolic invariant manifold S e = S for (13. ip for all e 
sufficiently close to 0, but the dynamics on S e may be different than on S. In particular, we 
are interested in the persistence of equilibrium points on S E near G. 
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4 Translational symmetry- breaking 

We now consider equations (13. ip in the case where e^O. Suppose ( c 



of (13. ip on the manifold S £ near 6, for < e <C 1. We use the local coordinates (I3.13P and 
write 



u(t; c, e) 



T a (t) 



Ei(a(t),c,e) 




where E = (E\, 0) T is a smooth function 

i?:lx [-c , c ] x [0, £ ] — > (Y 3 X R) n W 

satisfying 

sup{||£(a, c,e)\\yxr} < K < 00 V(c,e) G [-c , c ] x [0, e ] 
for some c > 0, e > 0. Substitution of (14. ip into (13. ip yields 



a \<pi + c 



e<(o) 








0(e) =c0x + L 



u(c) 



eL 



£1 




+ 



(4.i; 



T- a G(T a u*(0),0) 




7(u*(0)+cu;(0) + i;(c),c) 




^(«*(o),o) 





D u 7-(m*(0),0)m*(0) 







D u T(u*{0),0)v(c) 




/ V 



where 
q(a,c,e) 



-[J^w* (0) + cw*(0) +v + eE 1 ,c)- J"(u*(0) + cw*(0) + v, c)] - A^fV (0) + cu*(0) + v, c)£i 

+ j d u j ; -(«*(o) + c<(o) + w,c) j e; 1 - j d u j : -(«*(o),o)e 1 + 

+T- a G(T a (u*(0) + cu*(0) + Vl + eE 1 ),e) -T- a G(T a u*(0),0) 
is smooth and such that lim q(a, c, e) = 0. Projecting the above equation onto E using 

(c,e)->(0,0) 

the projection P and using (I3.14p leads to the following perturbation of the center manifold 
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equations (I3.12p 



where 



a = c — er(a) + e q(a,c,e) 
c = 



(4.2) 



and q(a, c, e) is smooth and such that lim q(a, c, e) = 0. 

(c,e)-K0,0) 

Remark 4.1 We make the following remarks concerning and ft4.3\ ): 

(1) IfQ were V equivariant, then it would follow that r (a) is a constant. Since it is assumed 
that Q breaks translational symmetry, then in general r is a smooth, non-constant 
function of a. 



(2) The function r(a) and the term q(a, c, e) in \4-ty are uniformly bounded in a for (c, e) 
near (0, 0). 

Proposition 4.2 For any e sufficiently close to 0, there exists a smooth curve O e which is 
the graph of a function c = e r(a) + e a (a, e) (where a is bounded and a (a, 0) = 0) such that 
a — in {4-2^ whenever (a, c) belongs to this curve. Consequently, Q e designates a curve of 
equilibria for the perturbed center manifold equations fiJJfy. 



Proof This is a simple application of the implicit function theorem. ■ 
From this proposition, we conclude that if c satisfies the relation 

inf {r(a) + cr(a, e)} < - < sup{r(a) + <r(a, e)} 

then a in ( 14. 2 p will become at at least one point o£l. Consequently, the flow line (a(t), c) 
for ( 14. 2 p contains at least one equilibrium point, which prevents the flow line from going from 
a = — ootoa = oo. The physical interpretation of this phenomenon is what is commonly 
referred to in the literature as propagation failure, or wave-blocking [H [121 [HI El [231 125] . 
On the other hand, if c satisfies the relation 

c c 

- > sup{r(a) + a (a, e)} or - < inf {r(a) + a (a, e)} 

e am ' ^ a6R 

then a in ( 14 .2p will be strictly positive, or strictly negative, but non-constant. Therefore the 
corresponding flow line (a(t),c) for (14. 2p will evolve at a non-constant drift speed, but the 
sign of this drift speed will not change on the flow line. Therefore, propagation failure (or 
wave-blocking) does not occur for this flow line. The situation is depicted in figure [3] 
We have thus proved the following 



17 



Theorem 4.3 Consider the abstract differential equation A3. 1\) satisfying Hypotheses \2.1\ 
IJ.il and \3.3\ For all e > close enough to 0, the semi-flow for this differential equation 
admits a two-dimensional normally hyperbolic C k+l -smooth invariant manifold S £ . Solutions 
of A3. 1\) which belong to So = S correspond to relative equilibria (with constant drift speed c). 
On the other hand, solutions of A3. 1\) which belong to S £ (for e ^ 0) correspond to perturbed 
relative equilibria, with non-constant drift speed. In the latter case, propagation failure (or 
wave-blocking) occurs if(c,e) belongs to a cone in parameter space, as illustrated in figure^ 



5 Parameter independent cases 

Hypothesis 13.3( a) implies that Z) c J 7 (m*(0), 0) ^ 0, otherwise L would also have (0, 1) T 
as an eigenvector. Thus, our results do not immediately apply to equations for which 
D c J 7 (m*(0), 0) = 0, which includes the cases where T does not depend explicitly on the 
parameter c, or on any other parameters. This comment is particularly relevant for the 
Fisher-KPP equation [T8] : 

du d 2 u _ . 

m = d^ + u{1 - u) 

which admits travelling wave front solutions U(x + ct) for all c in a neighbourhood of (£], 
even though the reaction term u(l — u) does not depend on any parameters. For cases such 
as this, the main ideas in this paper can still be used, but need to be slightly modified. 
First, the equation 

du 

— = Au + F(u) + eG(u,e) (5.1) 
at 

is studied without suspending the system (as was done in (13.11) ). We then assume hypotheses 
similar to Hypotheses l3.ll i.e. existence of a family of travelling wave solutions for (15.11) (when 
e = 0) parametrized by drift speed c. 
The linear operator L becomes 

L:Y — >Y 

L((p) =Acp + D u T(u*(0))cp 

with eigenvector 0i = £ u*(0) and generalized eigenvector 02 = u*(0). The rest of the analysis 
proceeds essentially the same as in section 3 and 4. In particular, one can show that there 



*For — 2 < c < 2, the wave front has an oscillatory approach to 0. Therefore, for many applications in 
biology, solutions with — 2 < c < 2 are not considered since they are not biologically relevant (the waveform 
attains negative values as it oscillates). 
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exists a normally hyperbolic two-dimensional center manifold S = {T a u*(c) ; a G K, c G 
(—c , c ) } C Y on which the flow of (15.11) at e = is given by (I3.12p . and which persists 
to a normally hyperbolic invariant manifold S £ = S for (15. ip when £ / 0. However, the 
perturbed equations on S £ become 

a = c — e ri(a) + e qi(a, c, e) 
c = er 2 {a) + e q 2 (a, c, e) 

where 

riWl = - ( /-Q)P ( r-MTM) j ^ = Qp ( W^") ) , (5 .3) 

and gi : 2(o,c, e) are smooth and such that lim 51,2(0, c, e) = 0. 

(c,e)->(0,0) 

For (c, e) near (0,0), the principal part of (15.21) is 

d = c — eri(a) 
c = er 2 (a). 

The dynamics of (15.41) can be quite complicated, depending on the functions r x (a) and 
r 2 (a). We will not give an exhaustive classification, but we will give a few illustrative 
examples. Many more examples are given in [23] for a specific three-species system of 
react ion- diffusion equations. In order to simplify things, we will suppose in our examples 
that the inhomogeneity is localized in the sense that both ri(a) and r 2 (a) vanish outside the 
interval I = (do, oi)- 

Example 5.1 Suppose the function r 2 (a) is of the same sign on I (without loss of generality, 
suppose r 2 (a) > 0). Then the phase space is separated into two regions by a curve C which 
is the union of two orbits of fl5.^[ ): one which tends to V = (a, c) = (do, £Ti(ao)) as t — > 00, 
and the other orbit which tends to this same point as t — > — 00, as illustrated in figure 
On one side of C, orbits with c ^ propagate from a = —00 to a = 00 (or vice-versa) 
without blocking, however the drift speed of the corresponding travelling wave is modulated 
by the inhomogeneity. On the other side of C, orbits of (5jj ) undergo a reflection by the 
inhomogeneity as they travel from a = 00 to the inhomogeneity, and then back to a = 00, 
as illustrated in figure The size of the interval of c-values for which orbits are reflected is 
approximately 0(e). 

Example 5.2 Suppose now that the function r 2 (a) changes sign over I . For the purposes of 
this example, we will suppose that r 2 (a) has two simple zeroes at a = p\ and a = p 2 in (a , ai), 
and suppose that r 2 (a) > on (ao,pi) and on (p 2 ,a 2 ), and r 2 (a) < on (p 1; p 2 ). Then 4\ ) 
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(5.2) 



(5.4) 



has hyperbolic equilibrium points at (pi,er 1 (p 1 )) (node) and at (p 2 , eri(p 2 )) (saddle). A 
possible phase diagram is shown in figure® The stable and unstable manifolds of the saddle 
point separate the phase space into three regions. In two of these regions, orbits propagate 
from a = — oo to a = oo (or vice-versa) without blocking, whereas in the other region, 
propagation failure occurs by one of two mechanisms: either reflection off the inhomogeneity, 
or blocking by the stable node. 

Example 5.3 In Example \5.S\ above, instead if having a stable node, the node could be un- 
stable and surrounded by a stable limit cycle. The situation is depicted in figure [7| This 
illustrates another mechanism for propagation failure: the travelling wave tends asymptoti- 
cally to a state characterized by a waveform which oscillates about a steady state. 

As mentioned above, these three example only serve to illustrate some of the possibilities 
for the perturbed dynamics in (15.41) . An exhaustive classification would be quite tedious, 
and certainly beyond the scope of this paper. 

6 Conclusion and Discussion 

We have presented a model-independent analysis of the phenomenon of wave-blocking, show- 
ing that under some general hypotheses, this phenomenon is universal in the context of forced 
translational symmetry-breaking in parametrized families of differential equations on Banach 
spaces. The approach is general enough in scope to be applicable to reaction- diffusion partial 
differential equations, integro-differential equations, retarded partial functional differential 
equations, and more. Our main theorem to that effect is Theorem 14.31 

We note that from a practical point of view, explicit computations of the various terms 
in the perturbed center manifold equations (14.21) and (I5.2p require explicit knowledge of the 
various terms which contribute to it, for example the function w*(c), and the eigenvectors 0i 
and 4>2 of the linearization L. 

We have chosen to work in a Banach space setting so that our results could apply to 
as many situations as possible. In certain applications (e.g. when studying pulses which 
decay to in reaction-diffusion partial differential equations), it is customary to work in the 
Hilbert space £ 2 (M, W 1 ) (e.g. [151 123]). In this case, one can use the adjoint operator L*, 
the Fredholm alternative, and orthogonal projections for the various projections which were 
introduced in Section 3. 

Finally, we have discussed in Section 5 how our approach can be modified to study 
parameter-independent equations (such as the Fisher-KPP equation) which admit families 
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of travelling waves parametrized by drift speed. We gave a few illustrative examples of how 
localized inhomogeneities can effect the dynamics of travelling waves in these cases. 
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Figure 1: Schematic representation of the decomposition of Hypothesis 13.3( a) as it relates to 
the center manifold S, the family of steady-states G, and the coordinates (a, c) of the center 
manifold equations f)3.12p . 
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Figure 3: Flow of (14. 2 p on the perturbed center manifold S e . All flow lines have constant 
value of c, but a is not constant along the flow lines. Flow lines between the dashed lines 
undergo wave-blocking (by meeting equilibrium points belonging to Q e ). If r max and r min 
designate respectively the supremum and the infimum of the function r(a) on R, then the 
width of c- values for which wave-blocking occurs is approximately e (r max — r min ). 
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Figure 5: Flow of (15. 4 p in the case described in Example 15.11 Orbits on one side of C do not 
undergo propagation failure, however their drift speed is modulated by the inhomogeneity. 
On the other side of C, orbits coming from a = oo in a region of width ps 0(e) undergo 
propagation failure: they get reflected by the inhomogeneity and bounce back to a = oo. 
The two wide lines at c = on either sides of the interval / correspond to unperturbed 
steady-states c = in (13.121) . 
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Figure 6: Flow of (15. 4 ft in the case described in Example 15.21 The stable and unstable 
manifolds of the saddle separate the phase space into 3 regions. In two of these regions, 
orbits do not undergo propagation failure, however their drift speed is modulated by the 
inhomogeneity. In the other region (of width 0(e)), orbits either reflect off the inhomogeneity, 
or get attracted into the stable node. 
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Figure 7: Flow of (15 .4p in the case described in Example 15.31 Similar to the previous case, the 
stable and unstable manifolds of the saddle separate the phase space into regions in which 
orbits propagate along the entire real line, and regions where propagation failure occurs. In 
this stable limit cycle attracts certain orbits coming from a = — oo. 
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